Single Qubit Rabi Oscillation Decohered by Many Two-Level Systems 
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Recent experiments on Josephson junction qubits have suggested the existence in the tunnel 
barrier of bistable two level fluctuators that are responsible for decoherence and 1 /f critical current 
noise. In this article we treat these two-level systems as fictitious spins and investigate their influence 
quantum mechanically with both analytical and numerical means. We find that the Rabi oscillations 
of the qubit exhibit multiple stages of decay. New approaches are established to characterize different 
decoherence times and to allow for easier feature extraction from experimental data. The Rabi 
oscillation of a qubit coupled to a spurious resonator is also studied, where we proposed an idea to 
explain the serious deterioration of the Rabi osillation amplitude. 
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I. INTRODUCTION 

Quantum computing with superconducting circuits has 
become one of the most anticipated approaches for the 
actual realization of quantum computation. However, 
like other solid state qubit candidates, while enjoying the 
advantage of scalability and the case of manufacturing, it 
suffers from short coherence time. This is mainly because 
the qubit states used in these devices are macroscopic 
quantum states, in contrast to the microscopic states 
used in NMR, ion traps, cavity QED, etc., which ob- 
viously have fewer intrinsic degrees of freedom that are 
open to perturbation from the external environment. 

Based on recent experiments, it has been shown that 
the primary noise that dephases a superconducting qubit 
has an 1/f like spectrum in a given low frequency range. 1 
Although the mechanism of this process is still unknown, 
it is speculated to be related to the critical current fluc- 
tuations in the Josephson junction due to its couplings to 
many two-level systems^. There are several possible ex- 
planations. The one we will focus on is that the trapped 
charge on the impurities in the tunnelling barrier of the 
junction block the junction area by Coulomb repulsion, 
thus modifying the critical current and the Josephson 
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energ; 

There have been several theoretical studies with similar 
models used to investigate the decoherence by evaluating 
the perturbation on the Rabi oscillations of the qubhpii 
In this article we analyze analytically and numerically 
the details of the Rabi oscillations for a qubit weakly 
coupled to an environment of two level systems. The 
time dependence of decoherence reveals different charac- 
teristics in the short and long time regimes, which can be 
used to analyze experimental data and extract informa- 
tion about the actual qubit energy fluctuation. Instead 
of trying to solve the problem exactly at once, in Sec. 
II, we started from an approximated model of a qubit 
dephased by many spins, which contains all the essential 
features where we can understand them in Sec. Ill, IV, 
V qualitatively and quantitatively by both analytical and 
numerical means. Later in Sec. VI we introduce other 
elements such as exchange interactions between the qubit 
and the spins whose energy is close to the qubit 's. This 



results in the relaxation effect and sometimes the com- 
plete destruction of coherence when they are exactly on 
resonance, as explained in Sec. VII. 



II. THEORETICAL MODEL 

The theoretical model of the system we consider is basi- 
cally a qubit coupling to many other spins while undergo- 
ing microwave driven Rabi oscillation. The Hamiltonian 
is 
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7~tqb + ^spins + T^int + T^Rabi , 
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(1) 



The a and r in the Hamiltonian are Pauli matrices. £1 
and u>k are the energy splittings of the qubit and the 
spins. Ak are the qubit-spin coupling constants. 2a is 
the Rabi frequency and O' is the microwave frequency. 

In the charge trapping picture, the u>k are given by 
2 A — where A is the gap energy of the supercon- 
ductor and tk is the trapping energy of the impurity in 
a Hubbard model 8 . When > 0, the trap is attrac- 
tive. The LOk here represent the energy differences of the 
processes where a Cooper pair breaks up and one of its 
electrons hops onto the impurity. Such attractive im- 
purities have been found, for instance, in the so called 
_D x -centers in semiconductors formed by substitutional 
dopants in GaAs and AlGaAs alloy aV . We assume the 
density of the traps is sufficiently low so that they don't 
interact with each other. Also the junction barrier is 
thought to be so thin that electrons only hop onto one 
single trap instead of a chain of traps. 

Ak means the amount of qubit energy being changed 
when the /cth trap is occupied or deoccupied. The 7iR a bi 
represents the Hamiltonian of the qubit under the driv- 
ing microwave at frequency f2' with the rotating wave 



approximation. 

For the interactions between the qubit and the spins we 
only include the terms where the spins perturb the energy 
splitting of the qubit. Here we omit the qubit flipping in 
the first part of this article because we want to investigate 
solely the dephasing caused by the many-spin system. 
Since the charge trapping processes only modulate the 
Josephson energy, this can be viewed, for example, as a 
model for a Cooper pair box at the charge degeneracy 
point. However, this model is not limited to only this 
picture. It is generally applicable to all qubit systems 
with dephasing sources approximated as quantum two- 
level systems. 
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III. ANALYTICAL SOLUTION 

The Hamiltonian in Eq. is exactly solvable when 
the number of spins is infinity^. First we apply an uni- 
tary transformation, 



.n't 



U) k t 



[/(*)= exp [i— — 



^ 2 
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(2) 



So that in this rotating frame the effective Hamiltonian 
is 



H = — =— ^ + ]T A k r z k o z + aa* 



Define B = J^k ^ T k- Since commute with the Hamil- 
tonian, they are all constants of motion and their values 
have their own thermal distributions. The qubit energy 
splitting is modified by its coupling to the spins and be- 
comes f2 + 2(B). Since this new energy splitting is the 
one that is measured in experiments, to drive Rabi oscil- 
lation, we let ft' = Q + 2(B). Thus Eq. J3J is simplified 
to 



U = (B- (B))a z + ha x = B'a z + aa x 



The evolution of the physical quantity a z can be calcu- 
lated, 



(a z (t)) = Tr[e iK V z e- iH Vo] 
The initial density matrix is 



(5) 
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FIG. 1: Time evolution of the qubit decohered by many 
spins undergoing Rabi oscillation by a coherent microwave 
source, (a z (t)}. The Rabi frequency 2a/2ir = 200MHz, and 
8Q is generated from a special case where all spin energy 
splittings are chosen for simplification to be ujk/2ir = 1GHz 
and VXfc -4fc/ 27r = 50MHz at the temperature T = 150mK. 
For Sl/2n = 10GHz, these parameters correspond to 8SI/Q < 
0.005. Notice that the amplitude of the oscillation envelope 
is already down to 50% at t ~ 20ns, but it reaches 25% only 
after t = 80ns, while an exponential decay should reach 25% 
around t ~ 40ns. 



Because the initial state is separable, the calculation 
of Eq. © ends up with calculating terms like (B' n ) = 
(3) Tr[B' n po]. Since in this approximated model t% are all 
independent random variables governed only by temper- 
ature, for a very large number of spins, the value of B' is 
given by a Gaussian distribution with 



(B') = 0, 
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The quantity Sfl > is a measure of fluctuation of the 
qubit energy splitting due to finite temperature. We can 
then obtain the exact solution in a closed form 



(4) (a z {€))=-! + 
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sin 2 Va 2 + B'H 
a 2 + B' 2 6 



( 9 ) 

Although the spins considered here are coupled to a ther- 
mal reservoir which gives rise to the Gaussian distribu- 
tion of the variable B' , the same formalism still applies 
when the Gaussian distribution in Eq. © is replaced 
by any arbitrary distribution for other mechanisms that 
govern the fluctuations of the two-level systems, and the 
precise definition of <5il will have to change with it. 



so that the qubit is in the | J,} state and the spins are in 
thermal equilibrium at temperature T — l/fcg/3. This is 
close to the real state of the system when the temperature 
T <C n. Note that the time evolution operator can be 
rewritten as 



e mt = cos 



yV + B'H 
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IV. DECOHERENCE 

Since this integral in Eq. is difficult to carry out, 
first we look at a special case. When T — * 0, 5£l — > 0. The 
Gaussian distribution function goes to a delta function, 
which means no thermal fluctuation. In this limit 



lim (<j z (t)) = 

T^0 
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FIG. 2: Fourier spectrum o z (uS) near the Rabi frequency 
uj — 2a. We can use the width of the peak to estimate the de- 
coherence time of the Rabi oscillation. SQ> and 5Q < are the 
larger and the smaller values of two random 5Q. For larger 
8Q we have a wider width and thus shorter decoherence time 
and vice versa. For frequency d > 2a + 25Q the tail of the 
peak is dominated by the Gaussian term, therefore we can ex- 
pect a Gaussian decay of (cr z (t)} when t <C h/80.. While the 
frequency id is very near 2a, the peak goes like (uj — 2a) _1//2 , 
which implies a very slow decay at the long time limit. 



which is the ideal Rabi oscillation, because every spin 
is frozen to its own ground state. But as soon as we 
turn up the temperature, when the value B is allowed to 
fluctuate, the oscillation now has a decay pattern(see in 
Fig. (Q). Notice that the oscillation lasts much longer 
than an ordinary exponential decay but has a rather dras- 
tic decrease of amplitude in the beginning. 
The Fourier transform of Eq. © is 
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for \uj\ > 2a. The meaning of Eq. is clearer if we 
define uj' 2 = uj 2 — Aa 2 and rewrite it as 



a z (oj') 



1 
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\/uj' 2 + 4a 2 uj' 
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(12) 



Now the width of the peak uj' = 0, i.e. lo = 2a, is 
controlled by the Gaussian function with width 26Q, as 
shown in Fig. J2J. When <5f2 — * 0, the spectrum be- 
comes a delta function at the Rabi frequency uj = 2a. At 
high frequency the Fourier spectrum is dominated by the 
Gaussian term, which means that (a z (t)} has a Gaussian 
decay in the short time limit. The singularity at the Rabi 
frequency uj = 2a implies a very slow decay in the long 
time limit. The result agrees with the that obtained for 
Eq. © in the limit when a ^> Sfl, the oscillation envelope 
of (a z (t)) is given by cr z (0)[l + (2tSfl 2 / a) 2 }~ ? , where the 
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FIG. 3: Temperature dependence of the peak width and the 
decoherence time, (a) is the direct plot of Eq. (11311 using the 
same parameters as those in Fig. Q. (b) is the decoherence 
time, which is the inverse of (a). Both the width and deco- 
herence time saturate when ksT > huik- 



evolution begins with a Gaussian (quadratic) damping 
then changes to a slow power law decay of ~ ^a/S&t 11 . 

Since the mapping from uj to uj' is not linear, the Gaus- 
sian width 28Q in oj' is approximately translated into a 
a) in uj. Thus we can define the 



width 2{\l a 
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decoherence time as 
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It shouldn't be a surprise that the decoherence time of a 
Rabi oscillation depends on the Rabi frequency. In the 
limit of small Rabi frequency a <§; Sft, 1/T<f, — * 2i5fi, 
which means the slower the Rabi oscillation, the more it 
is affected by dephasing. Also the physical meaning of 
is similar to its counter part in an exponential decay, 
T 2 * , since both of them are the product of pure dephasing 
and do not involve relaxation. 

The temperature dependence of the peak width and 
the decoherence time is plotted in Fig. ©. Notice the 
saturation of decoherence time is due to the saturation of 
50, thus the spins with \ujk\ = |2A — £fe| < ksT is exper- 
imentally unfavorable. If the electron trapping picture is 
correct, then better tunnel junction barrier materials or 
fabrication procedures should be sought out that doesn't 
contain electron traps that have attractive trap energy 
close to the superconductor gap energy. 



V. NUMERICAL SIMULATION 

To investigate more realistic situations we are required 
to solve the time-dependent Schrodinger equation numer- 
ically. Since the dimension of the Hilbert space grows 
exponentially with the number of spins, it is necessary 
to use an efficient and accurate method to propagate the 
dynamics of the system for long times. We adopt the 
Chebyshev expansion methodic to carry out this task. 
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FIG. 4: Numerical simulation for the Rabi oscillation of one 
qubit coupled with 14 spins at T = 200mK and T = lOmK. 
The left column are the real time Rabi oscillations, and the 
right column are their Fourier transforms. The dotted lines 
in the right column are the exact Fourier transform obtained 
in Eq. lilt for the limit of infinite number of spins. T^, cal- 
culated using Eq. 1131 . are 84ns and 165ns respectively. No- 
tice that the small bumps in the oscillation envelope in the 
T = lOmK are due to the effects of finite spins and some spins 
are frozen. 



FIG. 5: Rabi oscillations of a qubit dephased and relaxed by 
a many-spin system. The parameters are the same as those 
used in the T = lOmK graph in Fig. Q. The only additional 
parameter, the relaxation time Ti, calculated from Fermi's 
golden rule is 1/is. The solid lines are the numerical results, 
and the dashed lines are the approximations in Eq. 11511 and 
Eq. . 



VI. DEPHASING PLUS RELAXATION 



For initialization, a randomly generated state is propa- 
gated with the Hamiltonian Ti q b + TC S pin in an amount of 
imaginary time, i/3, instead of real time, and then normal- 
ized, so that the initial states have a Gibbs distribution. 
Later the normalized initial state is propagated in real 
time with Hamiltonian Ti., or equivalently Ti. 

First let us start with the simplified model previously 
considered analytically. The coupling parameters \A k \ 
are randomly chosen between 10 ~ 30MHz, and the \u)k\ 
of total 14 spins are uniformly distributed per decade 
from 2MHz to 10GHz. Since the chosen qubit bare en- 
ergy splitting is Q, = 2ir x 10GHz, much higher than the 
temperature, the initial state here is basically the same 
as Eq. ©. 

Fig. @ shows the results at T — 200mK and T = 
lOmK. The graphs in general look like Fig. Q and 
Fig. P|l. For the same set of parameters, higher temper- 
ature obviously gives faster and larger reduction of the 
oscillation amplitude. From the graph we can see that 
the right wing of the peak in the Fourier spectrum actu- 
ally consists of many small peaks due to the finite amount 
of spins. The physical meaning of T^, illustrated in the 
real time evolution graphs, is to characterize the process 
of the initial reduction of the amplitude, in this non- 
interacting spin case, a Gaussian decay. The much slower 
power-law decay that follows is just a general symptom of 
the coupling to many-spin systems and contains the infor- 
mation about dephasing only in its amplitude which was 
left over by the initial reduction. The rate of power-law 
decay cannot be used to extract any useful information. 
This special feature may well explain part of the reason 
for the problem of seriously reduced Rabi oscillation am- 
plitude which plagues many experiments. 



The relaxation of the qubit by a many-spin system can 
be described by adding a term to Ti- ul t , so that now 
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A k T z k o z 
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This kind of processes has been carefully studied by Pal- 
adino et. aZ.-. The relaxation effect caused by spins with 
energy splittings close to the qubit's alone in general gives 
the Rabi oscillation an exponential decay through the re- 
laxation induced dephasing. 

The simulation of both dephasing and relaxation 
present can be done easily in our program. The newly 
added spins remain non-interactive among themselves. 
Fig. (JSJ shows a simulation result with the relaxation 
time Ti S> T^^. In the real time evolution graph we can 
clearly see the three stages of the decay process, which 
starts with a fast Gaussian decay followed by slow de- 
cay, and then later the exponential decay finally takes 
over. Similar behavior has also been observed in the re- 
sult produced by a qubit under the direct influence of 
1/f noised. In the graph of the Fourier transform of the 
same data, the original sharp peak in the Rabi frequency 
is now smeared. Since this problem cannot be solved 
analytically, we made the approximation of multiplying 
the oscillating part of Eq. © with an exponential factor 
e -t/2Ti SQ foak it becomes 



(15) 
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The Fourier transform of this new function is 
f°° 1 

(16) 

and Re[<r z (w)] has been used in approximating the nu- 
merical result in Fig. ©. Although the approximation 
deviates from the numerical data significantly in the fre- 
quency domain, which is particularly obvious in the tails 
of the peak, they are in reasonable agreement in the time 
domain. This approach thus allows us to easily extract 
the relevant physical quantities such as and T\ from 
experimental data in a system of such complexity. 



where f2 res denotes the frequency of the two-level res- 
onator, which couples to the qubit with coupling constant 
g. The resonator is also dephased by the same group of 
spins that dephase the qubit. Here we drop the relax- 
ation effect, because it is less significant. As in the single 
qubit case, this model is exactly solvable in the limit of 
an infinite number of spins. After applying the unitary 
transformation 



/ n't , n't , M-^Lj k t A 

U = exp I i— a* + i—a* es + i ^ -j-r§ \ , (18) 
the Hamiltonian becomes 



VII. COUPLING TO SPURIOUS RESONATORS 

Recently experiment by Simmonds et. a& reported 
the observation of qubit coupling to some spurious res- 
onators. When the qubit is in resonance with one of these 
resonators, the Rabi oscillation pattern has a significant 
decrease in amplitude with irregular shape but not in de- 
cay time. In general, when a qubit couples with another 
two-level system while undergoing the coherently driven 
Rabi oscillation at its original frequency, it should dis- 
play beats in the patternLii instead of a decrease in the 
amplitude. As previously discussed, the reduction of am- 
plitude is due to dephasing, thus we suspect the reason 
for the amplitude reduction is that the spurious resonator 
itself is also being dephased, which further decreases the 
amplitude and smears the beats. 

The theoretical model we consider is as follows, 

7~L = W q b + Hres + Wspins + nint + T^Rabi ', 

n z 

rtqb = , 
nj _ ^rcs 2 

es _ 2 rcs ' 
_ Y^wfe 

T-spms — / , 2 k ' 
k 

Hint = g(a + a- cs + <j-<7+ s ) 

k k 

H Rabi = a(a+e- m,t + a-e in,t ), (17) 
I 



W = 2 e + — — a z cs + g{a + a Iea + a a+ s ) 

+ Y. A ^t° Z +J2CkT z k a* es + aa x . (19) 

k k 

Because it is the renormalized frequencies of the qubit 
and the spurious resonator that are in resonance at the 
driving microwave frequency n' , we can obtain an ef- 
fective Hamiltonian by defining B rcs = ^2 k CkT§ and 
B' Ies ee B res - (B les ), to obtain 

H = 9(<T + a- 3 + a-a+J + B'a* + B' Ies a* Ies + aa x . (20) 

This is a 4 x 4 matrix that can be easily diagonalizcd, 
though the result is a little messy. By the same token, 
physical quantities such as (cr z {t)} can be calculated by 
integrating out all possible values of B' and B' les . 

The density matrix of the initial state is chosen so that 
both the qubit and the spurious resonator begin with the 
| |) state, 

po = I | qb ® 1 U ® e-ez* V*. (21) 

The analytic result of the Fourier transform is approxi- 
mately 



~o z (lo) ~ dt dB' dB' Ies r(B',B' Ies )e e t[uJ±f(B ' > B '^ )] \ (22) 

J J — oo J — oo 

f(B',B' Ies ) = 2yfo? + B' 2 + B' 2 S + g 2 /2 ± v^/4 - 2B'B' Ics g 2 + AB' 2 B' 2 S + g 2 a 2 + AB? cs a\ (23) 

I 

where the definition of Sn rcs is similar to Eq. JHJ, which is the energy fluctuation of the spurious resonator. The 




FIG. 6: Effect on Rabi oscillations caused by a dephased 
spurious resonator. The solid line represents the case where 
the qubit is coupled to a spurious resonator, and the dotted 
line is not. The parameters for both are 2a/2n — 100MHz 
and 8Q./2-K = 145MHz. Those for the coupled resonator are 
g/2-K = 15MHz and 8Q tcs /2tt = 30MHz. Both qubit and 
spurious resonator couple to the same 14 spins. Notice that 
because <5f2 > h, the initial reduction of amplitude ends even 
before one period of Rabi oscillation. 



function r(B', B' les ) is the other part of the integrand that 
does not depend on <5f2, 5Q ICS , or t. The variables B' and 
B' Tes are in fact not independent as suggested in the way 
the integral is written, which is only for the purpose of 
easier explanation. 

In the case when <5f2 rcs <C <5f2, we can see that in the 
frequency exponent f(B', B' Tes ) in Eq. the first three 
terms in the square root are similar to the frequency 
exponent in Eq. besides an extra term B'^ cs . This 
means that the Rabi frequency fluctuation is slightly in- 
creased by the fact that the resonator is being dephased. 
The rest of the terms in the square root are causing the 
beats. If we freeze B' ICS = for these terms, they would 
not fluctuate, and the oscillation would display a clear 
beat pattern. But if B' [es is allowed to fluctuate, the 
beat frequency will also fluctuate, which smears the beat 
pattern. This suggests that the disappearance of beats 
is due to the dephasing of the spurious resonator itself 
when 5fl ICS > g. However, notice that the argument 
above applies only when the resonator fluctuates mildly, 
5fl lcs < SQ; otherwise the resonator would be brought 
out of resonance with the qubit. 

Fig. © shows the comparison of the numerical re- 
sult of a qubit dephased by 14 spins with and without 
a spurious resonator, which is also dephased by the same 
group of spins. We can see that the Rabi oscillation be- 
comes somewhat irregular with reduced amplitude while 
the long-time slow decay still persists without the sign of 
fading away. This explains the experimental observation 2 
that decoherence time is not reduced by the coupling to 
the spurious resonator is because the time has already 
passed and the Rabi oscillations have entered the fea- 
tureless slow decay regime. Notice that the irregularity 
starts to appear only after t ~ 2tt/<7, which suggests the 
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FIG. 7: Simulation of the spectral probing of resonance on the 
qubit-spurious resonator system. In both graphs, the vertical 
axes are the value of (<t z ) at the steady state, and the horizon- 
tal axes are the driving microwave frequencies. The spurious 
resonator here has an energy splitting at 10GHz. The left 
graph shows when the qubit energy is detuned from the res- 
onator at 9.95GHz, the most visible peak is the qubit and the 
resonator peak is barely visible. When the qubit energy is 
tuned close to the resonator frequency, as shown in the right 
graph, level repulsion takes place. Notice that the spectral 
peaks here are clearly seperated even though 8Q, and <5f2 rC s 
are both greater than g. 



cause is due to the smearing of the beats. 

However, with the parameters used in Fig. (JfJJ where 
<50 and SQ rcs are both greater than g, one might suspect 
that since the fluctuations of the energy level splittings 
are strong enough to smear the beat pattern in the Rabi 
oscillations, they could also destroy the avoided crossings 
observed in the spectroscopic data of the qubit transition 
frequency. Here this situation is numerically simulated 
and the result is shown in Fig. J7J, where the parameters 
are mostly the same as those in Fig. except a much 
weaker driving microwave power of 2a/2-n — 1MHz to 
ensure that the probability of excitation to the up state 
is linear to the driving power and additional relaxation 
spins for faster convergence. We show that the spectral 
peaks in the simulation are still clearly seperated and 
that the dephasing by many two-level systems is funda- 
mentally different from those by harmonic baths. 



VIII. CONCLUSION 

We have analytically and numerically solved the model 
of a qubit undergoing Rabi oscillation while being de- 
phased by a non-interacting many-spin system. It is 
found that the oscillation pattern, in the short time limit, 
has a Gaussian decay, because the spins have no inter- 
action among themselves. In the long time limit, it is 
fundamentally different from a qubit coupled to a heat 
bathi^, for its longer than exponential decay. 

The Fourier transform of (a z (t)) has a special shape 
and thus requires a definition of its spectral width. We 



7 



define the dephasing time in Eq. (|13(l to quantify the 
simple relation between dephasing time, the Rabi fre- 
quency, and the fluctuation of qubit energy splitting. We 
found that the dephasing by a many-spin system governs 
the initial amplitude reduction of the Rabi oscillation, 
which ends around t ~ T$. The slow decay later is a gen- 
eral behavior of a qubit coupled to a many-spin system. 

When the relaxation mechanism is added to the sys- 
tem, numerical simulation shows that one more stage of 
exponential decay appears after the slow decay in the 
pure dephasing case. Since this model cannot be solved 
exactly, we established an approximation method to al- 
low for easy extraction of physical quantities of decoher- 
ence times from experimental data. 

We have also numerically reproduced the serious reduc- 
tion of Rabi oscillation amplitude caused by a spurious 
resonator by considering the situation that the spurious 
resonator is also dephased by the same group of fictitious 



spins. The reduced amplitude of oscillation and the per- 
sistence of slow decay are both demonstrated. 
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